Chaos and residual correlations in pinned disordered systems 
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We study, using functional renormalization (FRG), two copies of an elastic system pinned by 
mutually correlated random potentials. Short scale decorrelation depend on a non trivial boundary 
layer regime with (possibly multiple) chaos exponents. Large scale mutual displacement correlation 



behave as \x 



the decorrelation exponent /i proportional to the difference between Flory 



(or mean field) and exact roughness exponent f. For short range disorder /i > but small, e.g. for 

random bond interfaces /i = 5£ — e, e = 4 — d, and fj, = e( 35 1) for the one component Bragg 

glass. Random field (i.e long range) disorder exhibits finite residual correlations (no chaos [i = 0) 
described by new FRG fixed points. Temperature and dynamic chaos (depinning) are discussed. 



Low energy configurations in glass phases induced by 
quenched disorder are non trivial and their statistics are 
far from fully characterized. It was proposed Q, 0, Q 
that they exhibit "static chaos" i.e. extreme sensitiv- 
ity to a small change of external parameters, denoted S 
(e.g. disorder, temperature, interactions). Suitably de- 
fined correlations, or overlaps, between two (sub)systems 
are argued to decay to zero beyond the overlap length 
L«5, which diverges at small S as L$ ~ (5 _1 / Q , a the chaos 
exponent. A phenomenological droplet description pre- 
dicts 0] a = 4f — 6 where —6 is the thermal eigenvalue 
and d s the droplet (fractal) dimension, which should ap- 
ply to spin glasses as well as to disordered elastic systems 
(there d s = d) 0j- Numerical simulations Q have qual- 
itatively at least confirmed this picture for systems of 
low dimensions d. In higher d and in mean field models, 
the presence of (temperature) chaos is far more elusive 
and it is still highly debated whether Lg is infinite, or 
extremely large Chaos is a central issue for the the- 
ory of glasses and can be explored in many experimental 
systems either by varying external parameters or con- 
structing actual correlated copies [E 0> HI O ne wants to 
develop analytical tools beyond mean field for this prob- 
lem, and also investigate residual correlations at large 
scale. 

Here we focus on the simplest class, elastic manifolds 
in random potentials. Upon an applied force / they ex- 
hibit a depinning transition, also with non trivial rough 
configurations 14] . They have numerous experimental 
applications for pinning of e.g. domain walls in magnets 
|9j , electronic crystals, density waves lfj , vortices in su- 
perconductors [111 1 1 2| , and wetting 13j . They exhibit 
energy dominated glass fixed points where temperature 
is irrelevant, Tl ~ L~ 6 , 9 = d — 2 + 2£ the free energy 
fluctuation exponent, £ the roughness exponent of the 
pinned configurations u ~ . Chaos in pinned mani- 
folds was studied mostly via scaling arguments . The 
directed polymer (d — 1) was studied numerically and 
via analytical arguments for N = 1 indicating a = 1/6 
in agreement with droplets and recently, on hier- 
archical lattices In d = 2 chaos was demonstrated 



for periodic systems near the glass transition T g using the 
Cardy Ostlund RG, and linked to universal susceptibility 
fluctuations as in mesoscopic disordered metals [l8| . 

Pinned manifolds can be studied using functional 
renormalization (FR G) within an expansion in e = 4 — d, 
to one loop [l^. [2(J and recently to higher orders [2l| . 
In the large N limit [13 it matches results from mean 
field theory |2jj. The main aim of this paper is to inves- 
tigate static chaos using the FRG. It also opens the way 
to investigate "dynamic chaos" in driven situations, e.g. 
at depinning / = f c , or in moving phases [24|. where no 
other analytical tool seems available. Results also pro- 
vide a test for the FRG, as a complete understanding of 
the field theory (FT) at low temperature remains a chal- 
lenge. In particular we unveil a chaos boundary layer 
structure of the FT encoding for switches between low 
lying states, which parallels the one found to encode for 
rare thermal fluctuations [2^ |2|| . 

We first investigate how two identical copies of a dis- 
ordered elastic system decorrelate when they experience 
slightly different disorder at T = 0. Consider the hamil- 
tonian: 



H = 



d a x[-(Vu 1 ) 2 + V i (u l (x),x)] (1) 



where the two copies i = 1,2 are not mutually inter- 
acting. They feel two different but mutually correlated 
random potentials taken, with no loss of generality to be 
gaussian with a correlation matrix: 



Vi(u,x)Vj{u',x') =Ri j (u-u')6 d (x-x') (2) 

and Rij(u) = Rij(-u). We denote the forces Fi(x,u) = 
d u Vi{u,x). This is thus a two layer version of the stan- 
dard pinned elastic system, and the only "coupling" is 
due to the statistical correlations. For a short range 
(SR) function R(u) it models an interface with ran- 
dom bond (RB) disorder, for a long range (LR) function 
R(u) ~ a\u\ at large u it describes an interface in a ran- 
dom field (RF) of local variance a. For a periodic R(u) 
it models a pinned one component CDW or Bragg glass, 
with generalizations to A^-component displacements u a . 
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We start with the case where all Rtj (it) behave similarly 
at large u, e.g. Vi — V ± WS, V and W being uncor- 
rected and identically distributed. The correlations are 
characterized by the matrix ((..) are thermal averages): 

Cij (x - x') = < (u*(x) - u l {x')){u3(x) - u3(x')) > 

The behaviour usually assumed from scaling arguments: 

C 12 (x)=x 2 <f(6x a ) (3) 

corresponds to an overlap length L$ ~ S" 1 ^". Here we 
show how to obtain from FRG and compute the decay 
at large distance: 

Ci2(a:) ~x-oo a: 2 *"" (4) 
We study the statics using the replicated Hamiltonian: 

ai abij 

One performs either Wilson RG varying the short scale 
momentum cutoff A; = Ae~', or, equivalently to one 
loop, compute the effective action of a uniform mode 
adding a small mass (i.e. a term \rv?u\ in with 
to/ = me~ l . Defining the rescaled dimensionless correla- 
tor through Rij(u) = A^mj Rij(um^), with 1 = 
87r 2 , one derives the T = one loop FRG equations: 

diRijiu) = (e - 4C)fly(«) + CuR'^u) 

+~^» 2 - ^{^{R'lM + R"m ( 6 ) 

We also use their dynamical version, in terms of random 
force correlators Ay- = —R'/j, which, to one loop, de- 
scribes also the depinning threshold f = f c , via standard 
dynamical techniques |2jJ, l25| . Considering identical dis- 
order in each layer, A(u) = Ai 2 (u) satisfies: 

diA(u) = (e - 2C)A(«) + C«A'(u) 

-A'(u) 2 - A"(«)(A(u) - A(0)) + TA"(u) (7) 

and we have denoted the positive quantity 

f = A u (0) - A(0) = l( Fl (0,x)-F 2 (0,xW (8) 

measuring the difference between the disorders in the two 
layers (at the bare level Ti—o ~ 5 2 ). An(a) = A 22 (u) 
obey identical uncoupled equations, same as (0 setting 
T = 0, and is known to develop a cusp and become 
non analytic beyond the Larkin length R c — e lc /A. T 
in (JSJ formally plays in the FRG equation a role rem- 
iniscent to a "temperature" in rounding the cusp at 
u = 0. Its flow under RG and its physics are how- 
ever very different from temperature and must be here 
determined self consistently from 1)718(1 . It is rather 



a static random force, which encodes static displace- 
ments below some scale, and broadens the interlayer 
correlator preventing it from becoming non analytic: 

A 12 K 4 - /a) -> hih-hf^iA ~ 4) with 

fi = Fi(0,q)/q 2 (in Fourier). This term arises because 
the statistical translational invariance (of disorder under 
u a( x ) ~ * u a( x ) + fi( x )) present in each layer [18|. does 
not hold for mutual displacements. There is no mutual 
pinning, the two layers are pinned independently, but 
there is some residual correlation |27j . 

Let us analyze the flow equation Q for very small S 2 ~ 
Tq. Below R c the non linear terms are unimportant: the 
unrescaled A(0) and A n (0) do not flow, f t = e^~ 2 ^% 
and A(u) and An(u) remain approximately equal. Very 
near R c , when TiA' n (0) ~ eAn(O), then A(u) starts to 
differ from An (u) within a boundary layer (BL) region of 
width u ~ T; c -C 1 around u = and remains analytic. 
Outside of this BL region, for u ~ O(l), A(it) is still 
almost equal to An (it) which then converges to the non- 
analytic FP function A u (u). The flow of A(u) beyond 
i? c , and its BL form, is much more difficult to study. For 
the case at hand (all R^j with similar large u behaviours), 
it is clear that TJ keeps growing slowly. Two important 
questions are thus (i) how to characterize this growth, i.e. 
the overlap length (ii) where does it lead to, i.e. whether 
there is a finite T fixed point for A(w), i.e. are there 
residual correlations. We first address (ii) and start with 
the periodic problem. 

1- Periodic problem (CDW, Bragg glass): one assumes 
that fixed points Ahj(ii) = A n (u) = A^u) = - 
\u(l — u) are reached within each layer. Using [2s|. one 
finds that there is no non-trivial T > fixed point. Thus, 
A(u) = A 12 (u) flows slowly to zero as: 

d t A(u) = eA{u) + (e/36)A"(u) (9) 

The leading decay, using potentiality J Q A(u) — 0, is 
Ai(u) — gi cos(2nu) with digi — e(l — )g;. The cor- 
relation @ results from the FRG flow of A from repul- 
sive to attractive FP and the general formula Cij(q) = 
Aij.i(0) I (q 2 + m 2 ) 2 , exact for q = 0, m = mi and to O(e) 
(one loop) accuracy for q ^ m (q — Ki, m ^ Q for Wilson 
RG). At large scale = for the periodic problem): 

C 12 (q) « A^O^-V^-^l^i) ~ g -W (10) 

Thus the correlations of the displacements in the two 
samples decay as in Eq. (J3J with: 

^ = £((<Lt 2 /36) - 1) + 0(e 2 ) (11) 

For general periodic problem, e.g. the higher component 
Bragg glass, to 0(e) one finds (i = e(A* af3 (0)K°K% - 1) 
thus n = rjK ~ e i where r\K a is the Bragg glass expo- 
nent for decay of translational order in a single copy 
< e %K ° u *e~ iK ° u x ' > ~ x~ VK o . One can compare the 
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prediction (|llfl with the one from Ref. 18] in d = 2 
near T — T g . There, disorder generates an additional 
term SH rep = =± J x [crVi^Vuf + 2aVul Vit^] in J5J), with 
o\ growing as Bl and <j\ — > a*. This implies that 
C u {x) ~ Bln 2 |x| while C u (x) ~ a* In | x\, i.e. no ab- 
solute decay of interlayer correlations (/i = 0), although 
decay relative to intralayer ones. A numerical study |28j 
is in progress to determine if a* > at T — |2flj . 

Short range (random bond) disorder. There also we 
find that A(it) flows to small values. Thus we study the 
linearized part of H7I8|I around A(it) — 0. It is equivalent, 
and simpler, to study directly the unrescaled linearized 
RG equation for R = R 12 , diR(u) = -e d R'( l {Q)R"{u). 
Since each layer reaches a fixed point, i.e. Ru(u) — > 
Rti(u), one has ^(0) - e^ e + 2 ^ 1 Rl'{(0). Thus fl,(u) 
converges to a diffusion front (i?; = o(w) being peaked): 

~ -™=e-™ /(4Di) / dit'W) (12) 

with A = jJtfl'eW - e 2 ^. Thus A,(0) - -i*f(0) - 
D;" 3/2 - e~ 3C ' and from (JTUJl one finds (£j with: 

M = M^s = 5C - e (13) 

always strictly positive (to one loop Crb = 0.28e), i.e. 
the FP A = is attractive. This is easily generalized 
to the 0(N) model with the result n SR = (N + 4)C - e, 
replacing R"i(0) by d a dpRu(0) ~ 5 a p and considering 
diffusion in a TV-dimensional space. Finally, one shows 
that connected correlations of the free energies F, of each 
layer are given by F\F) C = R l3 {0)L d 30]. From JT2J) we 
find that the decay of energy correlations is also of the 
form F\F^ ~ L 29 -^ B . 

It is noteworthy that the replica variational method 
|23l ] yields Flory exponents, i.e. fj, = 0. The FRG cap- 
tures the physics of the decorrelations beyond the Flory 
approximation. Since (|13|) results from linear FRG at 
A* = 0, one may look for simple interpretation. Since 
each ground state satisfies V 2 u l (x) = Fi(u l (x), x), resid- 
ual correlations in position can be interpreted from cor- 
relation in pinning force resulting from accidental en- 
counters u 1 — u 2 (i.e. averaging V 2 ^, 1 ^ 72 ^ 2 , neglecting 
all correlations on the r.h.s. and using u 1 ~ yields 
i|4H3|0 . This argument (see also 0), is delicate how- 
ever, since dimensional estimate (Flory) fail here, as well 
as attempts to make it self-consistent. In fact, linear 
analysis still involves non trivial A l:L fixed point and ljl3|l 
is confirmed by FRG [3lj |. 

3- Long range (random field) disorder. A different sce- 
nario occurs there. Using now the linearized unrescaled 
A(u) equation, i.e. <9;A(w) = e el A n (0)A"(u), yields 
a diffusion form for A;(u) with now Di ~ e 2 ^', thus 
fi = 3(— e. Since for the RF case ( = e/3 to all orders (i.e. 
Flory is exact), the flow near A = is marginal along a 
line of FRG fixed points. These new FP describe the 



residual (i.e. large scale) mutual correlations, which we 
now describe. In the RF case, remarkably, the fixed point 
equation can be integrated for any fixed T . Also the 
values of / + °° duAij(u) = —RL(u = +00) = ov, are in- 
variants (exactly to all orders) of the FRG flow, where 
we denote an — a and a 12 = cr' = cr(l — S 2 ) the bare 
variances of inter and intra layer random fields. Defining 
rescaled variances A^o — e<r, Ada' = ea', the function 
7(7") = f dy^Jy — 1 — (1 + t) lny, a reduced tempera- 
ture r = 3T/(e£ 2 ) and the two lengths $ = and 

£ 3 = , the fixed point functions are then: 

V27(t)' 

A* 11 (u)= € -e Q yo(u/^) , A» = |eV(«/0 

y T {x)-l-{l + T)lny T (x)=x 2 /2 (14) 

and y T (Q) = 1. While A* u (u) is the standard RF T=0 
FP with a cusp non-analyticity at u = 0, A* (it) has the 
characteristic shape of a finite temperature FP rounded 
by a (Chaos) BL of width it ~ f. To find this FP one 
needs to determine T (i.e. t) for given RF variances 
a, a'. It is given by the self consistency condition JSJ, or 
equivalently £ 2 (1 + r) = £q. This yields: 

( r7^=(l + r)- 3/2 7 (r)/7(0) (15) 

which determines r given a' /a, from which we obtain 
£, then T and A* (it). One sees from (I15f) that there is 
indeed a fixed point A* (it) for each value of er'/er. In 
a sense the overlap length is infinite, since there is a fi- 
nite residual correlation (although there is a flow from 
infinitesimal T towards its FP value). It is described by 
universal functions, e.g.: 

Ci2{x)/C ll {x)^ x ^ 00 F d [a'/a\ (16) 

with Fd = (1 + t) -1 / 2 near d — 4, r being solution of 
(|15fl . It exhibits, to one loop accuracy, a small <5 sin- 
gularity, from (fT5|l , t ~ 8/ yJ\n(l/5). An interesting 
quantity to compute numerically is the ratio of center of 
mass fluctuations (adding a mass to, i.e. harmonic well) 
C u (q = 0)L£ji(<1 = 0) ^ m ^o G d [a'/a] with G d = F d 
near d = 4 32]. 

To discuss smaller scales (~ L§) we recall the droplet 
picture. Let u\ be the ground state (GS) of V\. A small 
SR perturbation V-x — V\ ~ 6 yields a new GS it2- For 

small system size L <C L$ u\(x) ~ U2(x) and the change 

1/2 

in GS disorder energy is ~ 

where Ei = J d d xV(x,u l (x)). The (small) probability 
p that another low energy state becomes more favor- 
able inducing a (rare) large switch u-i — U\ ~ L& is thus 
p ~ e /L e - SL<- d ^- e - 5L^r- ', as in thermal droplets. 
p becomes of order one at scale L$. How does the FRG 
capture this physics ? Consider the growth at short scale 
of (Fi - F 2 ) 2 - ARL 28 with AR = R n (0) ~ Ri 2 (0). 
From © one finds diAR = (d-26)AR- if 2 . This sug- 
gests that AR grows as 6 2 L d ~ unaffected by the non 
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linear term, until it reaches a FP value O(l), and repro- 
duces the droplet result at least for energy correlations. 

However, more complex things could happen, as 
seen on displacement correlations. The GS difference 
(u\ — « 2 ) 2 ~ TL 2t > is determined by the initial growth of 
T ~ L a . Within droplets T is proportional to p the prob- 
ability of rare GS switch, and a = adrop — d/2 — 9. Rare 
GS switches induced by thermal fluctuations was shown 
to be encoded in the Thermal BL |26|. Here, we expect 
that the Chaos BL (CBL) for A(u) describes switches 
induced by disorder perturbations (higher moments of 
responses being encoded in derivatives A( 2 ™)(0)). Wc 
found evidence for CBL solutions of the system (|7I8|I of 
the form A(u) = A(0)+f(f>{u/f)+O{f 2 ) and indications 
of a non trivial a < adrop- However, firm analytical as 
well as numerical determination of a is very delicate |33j . 
A non trivial a ^ adrop arising from the CBL would re- 
quire better understanding of minimization on all scales. 

That lj 718(1 yields a rich variety of behaviours is clear 
from e.g. the case of "subdominant chaos" Vi(u) = 
V(u)±v(u) where V(u) is e.g. RF disorder and v(u) RB 
disorder. Then the ground state should not change at 
large scale since 5L d / 2 < L 0RF with 6 RF = (d+ 2)/3. f 
should spontaneously flow to zero (correlations increase 
with scale) with adrop = d/2 — Orf = — e/6 and a (self- 
organized) CBL solution should form, similar to the TBL. 
Demonstrating such RG flow, and testing universal rela- 
tions with thermal moments in d = is in progress. 

Temperature chaos, i.e. two copies with different tem- 
peratures Ti, is described by Eq. iJBJ adding the term 
i(T- + f^R'l^ with fi ~ T,e- dl . The mutual correlator 
A still satisfies © with: f -> ±(Ti + f 2 ) + f where now 
f= i(A 11 (0) + A 22 (0))-A 12 (0). Even for identical bare 
disorder, one finds that although <9/T|; = o = one has 
dff |j =0 = -\(T 1 -T 2 ) 2 A'['{(Q) thus an effective non zero 
Tbare, and an independent short range An (it) — A i2 (ii) 
is generated at small scales from entropic difference be- 
tween the two copies. Apart from there being an addi- 
tional scale where the decreasing Tj cross the growing T, 
the flow of A12 at larger scales, and around L$ is similar 
(in the SR case) to disorder chaos. 

For dynamic chaos at depinning our result yields |3l| 
M-dep = 3£(2ep — e, with /idep > from two loop correc- 
tions [2lJ , and thus slow decay of mutual correlations of 
critical configurations at f c . However, it hold strictly 
in the fixed (or bounded) average center of mass ensem- 
ble, and a better treatment of the mode q — may be 
necessary (HH^. Eq. (|71gj) without the term A'(m) 2 
and e = £ = 0, describes the moving Bragg glass in 
d = 3. One finds (3(| residual correlations as in (JTSJ, 
with F d=3 « duA 12 (u) 2 /(2A n {0) 2 ) + 0{6 3 ) in terms 
of the bare disorder. 

To conclude we applied FRG to chaos in pinned sys- 
tems and obtained the residual correlations. In long 
range disorder these do not decay and are described by 



new FRG fixed points. The absence of chaos there is 
reminiscent of weaker or absent chaos in infinite range 
models of more complex glasses. Depinning and some 
quantum systems can be studied by this method. The 
self-organized chaos BL in the field theory uncovered here 
remains to be fully understood. We hope to stimulate 
numerical and experimental studies of correlations |3?| . 
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